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ABSTRACT 


A simple matching technique is explained which allows the computation of 
the response of the solar surface to perturbations which occur at any depth 
within the convective envelope of the Sun. This technique was applied to a 
perturbation of the convective efficiency (a-mechanlsm) , and of the non-gas 
component of the pressure (B-mechanlsm) In different regions of the convection 
zone. The results Indicate that either perturbation affects the solar lumin- 
osity. However, the a-mechanlsm has little effect In the solar radius, 
regardless of the location of the perturbed region, whereas the B~mechanlsm 
produces radius changes that become quite large If the location of the per- 
turbed region Is deep within the solar convection zone. 

Changing the mixing length ratio used In computing the solar convective 
envelope Is one mechanism capable of producing a fluctuation of the solar 
luminosity (refs. 1,2). However, there are other possible mechanisms which 
can also produce solar variations. To facilitate the Investigation of the 
effects of as many possible mechanisms as possible In an economical way, we 
shall describe here a simple matching technique for quick estimation of these 
effects. This technique cannot give the time evolution of the fluctuations, 
and the action of the perturbing mechanism must be confined to be above the 
bottom of the convection zone. However, this technique can give a rather 
accurate estimate of the amplitudes of the solar variations with very little 
computational effort. We shall first describe the mechanics of the technique, 
then Its Justification, and finally some applications. 

The procedure to follow Is very simple. First, one calculates the 
envelope structure of an unperturbed solar model using some standard stellar 
envelope code (cf. ref. 3). Then the perturbing mechanism must be Implemented 
In the evelope code, and a new envelope matched to the old envelope at a mass 
I point near the base, but within, the convection zone. By "matching", we mean 

I that the luminosity and the outer radius of the new envelope are changed until 

its pressure, temperature, and radius at the matching mass point become the 
same as those of the old envelope. After the fitting parameters are deter- 
mined, the differences between the new and old envelopes can be Interpreted to 
be the variation of the solar luminosity and the spontaneous change In the 
solar radius. 
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To justify this proceedurs, «ie need to answer the following questions; 

1. There are three variables: pressure, temperature, and radius, to be 

matched but there are only two fitting parameters - luminosity, and outer 
radius. Is It possible to match all three variables consistently? 

2. The matching point at which the matching Is done Is not a priori 
fixed. Are the results of matching Insensitive to this point? 

3. The envelope Integrations need to assume the luminosity to be uniform 
above the matching point. Is this procedure valid? 

4. The variation of the pressure, temperature, and radius at the 
matching point are taken to be zero. Is this a valid approximation? For 
how long? 

To answer questions 1 and 2, let us examine Table 1. This table shows 
the results of perturbing the solar envelope by Increasing the mixing length 
ratio o by -•■O.IZ from an Initial value of l.S. The first row at a given mass 
point (M^) gives the results obtained by matching the pressure and the 
temperature at M^., and the second row gives the results obtained by matching 
the radius and the temperature at As one can see, the variations In 

luminosity are Identical, and the variations In the solar radius all agree to 
within ISZ. The ratios between d log R and 6 log L are all about 6 x 10~^, In 
very good agreement with time-dependent calculations (refs. 4 and 5). 
Therefore, matching either pressure-temperature or radius-temperature give 
consistent results. Furthermore, this table shows that the results are 
Insensitive to the exact location of the matching point. 


For the luminosity L to be uniform In a perturbed envelope. It Is 
necessary that the time scale for the relaxation of a non-uniform distribution 
of the luminosity be very short. To estimate this time scale, we derived a 
diffusive type equation for a non-uniform distribution of luminosity AL In the 
almost-adlabatlc convective region (using mixing-length theory of convection) 
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where e - 4n pC Tr H, t ■ (4/3 a) (H/V). H Is the pressure scale height, V Is 
the convective ^loclty, and all other symbols have their standard meaning. 

The important time scales here is t which Is close to tlM eddy-tum-over time 
scale of the convective turbulence. Near the bottom of the convection zone, 
this time scale Is of the order of half a month. Therefore, convection can 
smooth the luminosity distribution In very short time. This smoothing time 
scale constitutes the lower limiting time scale for the applicability of our 
approach. 


The upper limiting time scale for the validity of this approach can be 
estimated as follows. An extra energy leaving the solar surface 5LAt Is 
approximately given by half of the gravitational energy released by a collapse 
of a layer of mass AH by an amount 6r in jradlus 
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This provides an estimate for the growth of the collapse 6r In time. 
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want the estimate of the variation In the outer radius 6R to be not complately 
washed away by the collapse fir, then we need fir to be leas than fiR* This 
Implies a restriction on the time period given by the Inequality 

At < 2 X 10^ W <yra) (3) 

where W • 6 log R/fi log L. In the case for perturbing a, W « 6 x 10"^; 
therefore At must be less than 120 yrs> 

We have answered all the previous questions on the validity of the 
approach. Now, let us briefly describe some results of applying the present 
technique to a calculation of the variations Induced by a perturbation of the 
magnetic pressure In the upper convection zone. Let 3 be defined as Pg|/Rt' 
where Is the material pressure, and is the total pressure, which In 

addition to Includes the radiation pressure and some "effective” magnetic 
pressure. Treated as a fluid, the magnetic field actually behaves quite 
similar to the radiation In the sense that Its polytropic Index Is also 3. A 
lowering of 6 would mean that the fraction of electromagnetic pressure Is 
raised relative to the material pressure. Figure 1 shows the results obtained 
by suddenly lowering 6 by O.IZ In the region between the solar surface and 
(the horizontal axis). The solid curve shows ^R/R^ and the broken curve 
shows fiL/L as a function of the depth of the perturbation (the matching point 
for all cases Is at * 0.99). One can see that as the perturbation goes 
deeper, the effects become more prominent. Especially, the perturbation In 
radius Increases substantially when 5B occurs below Mj, - 1-10”" the 
fi logR/fi log L ratio becomes close to 0.1. However, ve must stress that the 
present calculation Is very preliminary, and It can only give an order of 
magnitude guide to the general response of Che convective envelope to B 
perturbations. 
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TABLE 1. RESULTS FOR PERTURBING BY 


Ao + 0.1% FROM u - 1.5 


Matching 

Point 



6 log (L) 

6 log (R) 

A log (R) 
6 log {l) 

.9850 

.326(-3) 

.189(-6) 

.579(-3) 


.326(-3) 

.193(-6) 

.59K-3) 

.9900 

.326(-3) 

.189(-6) 

.580(-3) 


.326(-3) 

.185(-6) 

.566(-3) 

.9990 

.326(-3) 

.182(-6) 

.557(-3) 


.326(-3) 

.170(“6) 

.52K-3) 

.9999 

.326(-3) 

.185(-6) 

.567(-3) 


.326(-3) 

.165(-6) 

.507(-3) 






